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Experimental evidences are presented showing unusually large and highly anisotropic vibrations in the 
"simple cubic" (SC) unit cell adopted by calcium over a broad pressure ranging from 30-90 GPa and at 
temperature as low as 40 K. X-ray diffraction patterns show a preferential broadening of the (110) Bragg 
reflection indicating that the atomic displacements are not isotropic but restricted to the [110] plane. The 
unusual observation can be rationalized invoking a simple chemical perspective. As the result of 
pressure-induced d transition, Ca atoms situated in the octahedral environment of the simple cubic 
structure are subjected to Jahn-Teller distortions. First-principles molecular dynamics calculations confirm 
this suggestion and show that the distortion is of dynamical nature as the cubic unit cell undergoes large 
amplitude tetragonal fluctuations. The present results show that, even under extreme compression, the 
atomic configuration is highly fluxional as it constantly changes. 

High- density polymorphs of calcium show structural and electronic properties dissimilar to those of other 
elemental solids. Under compression at room temperature, Ca transforms successively, from the face- 
centered cubic structure (FCC, phase I) to body-centered cubic structure (BCC, phase II) at 19 GPa 1 ' 2 and, 
with further compression, to a simple cubic structure (SC, phase III) at 32 GPa 3 . Three more phases, beyond that 
of the SC phase stability field, were recently characterized from 100 to 300 GPa 4 . Unlike alkaline Group II Sr and 
Ba 5 ' 6 , an incommensurate phase in Ca was only found very recently at a pressure in excess of 200 GPa 4 ' 7 . 
Moreover, the SC structure, a rare occurrence in metallic elements, was found to be remarkably stable over a 
large pressure range from 32 to 100 GPa 6 . Density functional calculations employing generalized gradient 
approximation (GGA) have been performed in order to understand the stability of the SC phase 8 " 13 . It was 
concluded that the SC structure is dynamically unstable 8 and, instead, a tetragonal 1^1 amd structure 9 " 12 was 
found to be substantially lower in energy. More elaborate meta-dynamics 9 and molecular dynamics (MD) 11 
calculations suggested the SC phase could be thermally stabilized, although the technical validity of the calcula- 
tions has been raised 10 . Nevertheless, experiments by Gu et al 14 showed that at room temperature, even after 
sample annealing, the powder x-ray diffraction patterns could be indexed to an SC structure, irrevocably. In 
contrast, low temperature x-ray diffraction experiments have led to the conclusion of the existence of an ortho- 
rhombic Cmmm structure near 40 GPa at 4K 215 . Density functional GGA calculations showed that the latter 
structure is dynamically unstable and presents a higher enthalpy than that of the I^/amd structure 13 . Diffusion 
Quantum Monte Carlo calculations at GGA-optimized geometries, however, suggested that the relative stability 
of I^/amd and Cmmm phases at 50 GPa may be reversed but at low temperature the SC phase would still have a 
higher enthalpy 13 . A recent theoretical study 16 indicated that anharmonic vibrations might help to stabilize the SC 
structure 16 . However, in spite of extensive theoretical and experiment efforts, no satisfactory explanation on the 
electronic origin of the very flat and anharmonic potential energy surface and, equally significant, the reason why 
the enthalpy most stable I4i/amd phase has not observed are forthcoming. The discrepancy between theory and 
experiment is resolved and presented in the present work. Present results on low- temperature high-pressure x-ray 
diffraction experiments and first principles MD calculations suggest the s —> d hybridization of Ca leads to a 
dynamical lattice distortion of the SC structure rooted in an electronic instability. 

To resolve this problem, the pressure-temperature phase diagram of Ca was determined near the bcc —> sc 
transition. Powder x-ray diffraction experiments were performed using He as the pressure transmitting medium. 
Pure Ca samples (Sigma-Aldrich) were loaded in an inert gas atmosphere. Experiments were carried out at 
beamline BL10XU at SPring-8 using 0.4137A radiation, focused at the sample by a compound refractive lens 17 . 
Quasi-hydrostatic pressures were obtained from the calibrated spectral shifts of Al 2 0 3 :Cr 3+ luminescence and 
first-order Raman line of the diamond (anvil), measured and corrected at each temperature. Pressure and 
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temperature uncertainties were 0.1 GPa and 1 K, respectively. 
Several pressure-temperature paths were followed to explore system- 
atically the phase diagram with pressures ranging from 30-50 GPa 
and while varying the temperature from ambient to 14 K. 

Results 

The FCC — > BCC transition was observed between 19.3 and 
20.4 GPa at room temperature, in agreement with the reported 
value 1 . Then the Ca sample was cooled to 14K and further com- 
pressed up to 53 GPa. The onset of a structural transition of the 
BCC phase occurs at 39 GPa, which is much higher than the 
32 GPa observed under ambient temperature 2 . In a narrow pressure 
range above 40 GPa and at 14 K, a series of structural transforma- 
tions indicated by changes in x-ray diffraction patterns were 
observed suggesting the occurrence of possibly more than one inter- 
mediate phases (Fig. 1). At pressures above 42 GPa, the x-ray dif- 
fraction patterns became simpler (Fig. 1) and indexed to a SC unit 
cell. Hence at low temperature, Ca does not transform directly from 
the BCC to SC structure, but goes through a sequence of intermediate 
structural configurations. A distinct feature of the x-ray diffraction 
patterns is the full width at half maximum (FWHM) of the (110) 
Bragg reflection which is significantly larger than both the (100) and 
(200) reflections (Fig.l). The observed line broadening cannot be 
attributed to non-hydrostatic conditions, as the width of diffraction 
lines corresponding to similar d-spacings, measured at 14 K and 
300 K and around 40 GPa, differ marginally, indicating that He is 
still a good quasi-hydrostatic pressure transmitting medium. Also 
there is no evidence that the unusually large width of the (110) line is 
due to a splitting of the Bragg peak resulting from a symmetry 
change. A sequence of x-ray diffraction patterns was also recorded 



by warming the sample from 14 to 100 K, at around 42 GPa. No 
strong temperature dependence on the FWHM of the (100), (110), 
and (200) reflections was observed. In fact, the (110) line width 
reduced slightly from 0.67° to 0.64° (FWHM) with an increase of 
temperature from 14 to 100 K at 42 GPa, as measured to ± 0.01° in 
our experiments. In comparison, the FWHM for the (100) and (200) 
reflections remained to be 0.10° and 0.18°, respectively. Normally, 
the width of a Bragg reflection increases with the diffracted angle 18 . 
However, the exceptionally large line width and the fact that the 
(110) Bragg peak diffracts at a lower angle than the (200) reflection, 
rule out this explanation. The anomalous preferential broadening in 
the (hhO) Bragg reflections corroborates the fact that the Ca atoms in 
the [110] lattice plane are indeed disordered. Other possible effects 
like deviatoric stresses are not ruled out but considered unlikely given 
the peculiar line broadening has been observed repeatedly for pres- 
sures below 43 GPa and above 39 GPa at very low temperature 
(14 K) and not in other phases of Ca, namely phases II and III, which 
exist under close temperature and pressure conditions. Another 
usual observation is the line width of the (110) Bragg peak apparently 
decreases slightly with increasing temperature (Fig. 1). This anomaly 
will be discussed below (vide infra). 

The simple cubic unit cell parameters as a function of temperature 
along the 42 GPa isobar are calculated using the measured (100), 
(110), and (200) reflections. The results presented in Fig. 1 show a 
negative thermal expansion from 40 to 100 K. A negative thermal 
expansion is indicative of large anharmonicity in lattice vibrations 19 . 
Moreover, at temperatures above 80 K and 42 GPa, the estimated 
coefficient of linear thermal expansion of Ca of 2.1X10" 4 /K is an 
order of magnitude larger than that of the FCC phase at ambient 
pressure (2.23X10- 5 /K) 20 . 
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Figure 1 | (a) Selected x-ray diffraction patterns of Ca recorded at specified (P, T) conditions, (b) The measured line width for the (100), (110), and (210) 
Bragg reflections of the SC phase, (c) Variation of the unit cell with temperature at 42 GPa. 
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The peculiar behaviour of Ca at low temperature and high pres- 
sure gives rise to a different and confined field in the phase diagram, 
nested between that of phases II and III. This new stability field in the 
phase diagram is identified as that of phase III'. A schematic diagram 
for Ca is presented in Fig. 2. Within the phase III' field, at the lowest 
temperature reached in our experiment (14 K), a sequence of 
changes in the x-ray diffraction patterns was observed. None of the 
x-ray diffraction patterns could be successfully fitted to the predicted 
I4i/amd structure, a solution which has been definitely ruled out. The 
Cmmm unit cell identified earlier 215 as the best description of the Ca 
at 44 GPa and 4 K (close to the phase III-III' boundary) could not be 
confirmed as an unique and non- ambiguous structural solution to 
our x-ray diffraction results. It is plausible that small residual stress 
could help to stabilize the intermediate structures and the discrep- 
ancy from previous work could be due to the quasi-hydrostatic 
conditions. 

Discussion 

Why the structural distortion only occur in the [110] plane of the SC 
structure? This peculiar phenomenon can be explained from a simple 
crystal field picture. When compressed, alkali metals often transform 
between simple close packed structures 5 . At higher density, the near- 
est-neighbour atoms form confining potential localizing empty orbi- 
tals of higher angular momentum and then mixed with the valence 
state 21 . The s — > <i transition in Cs, first recognized by Sternheimer 22 , 
provided a good example of hybridization. In alkali and alkaline 
earth metals with low-lying d-orbitals the full realization of the 
s — > d transition is usually preceded by an intermediate pressure 
regime whereby unusually complex commensurate or incommen- 
surate structures are observed 23 . A common feature of these high- 
density structures is their tendency to adopt planar packing of 
atoms 5 ' 21,23 ' 24 . The participation of d-orbitals results in the local- 
ization of electron density in the interstitial region. This point of 
view, although empirical, has helped to rationalize the structural 
features, transformation sequences and charge distribution in Rb 
and Cs 21 . 

A previous calculation indicated that the d- character of the val- 
ence band increases substantially accompanying the BCC — > SC 
structural transformation in Ca 11 . Therefore, it is reasonable to 
assume that the promotion of the 4s electron into the 3d orbital of 
Ca is the cause for the transformation. Since each Ca atom in the SC 
unit cell is surrounded by an octahedral arrangement of nearest 
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Figure 2 | The schematic phase diagram of high pressure Ca. Phase 
boundaries are the results of the many temperature and pressure paths 
explored in this study. 



neighbours, under the crystal field 25 of the surrounding atoms, the 
d-orbital manifold splits into triply degenerate t 2g and doubly degen- 
erate e g sets of orbitals. As there are at most two valence 4s electrons, 
the promotion of an electron from 4s to 3d (4s 2 — » 4s3d) gives rise to a 
2 D state. Double excitation (4s 2 — » 3d 2 ) will result in an 3 F state. Both 
electronic states are orbitally degenerated and, according to the Jahn- 
Teller (JT) theorem 26 , a lower- energy structure can result from the 
distortion of the local octahedral environment to a lower symmetry. 
The lattice distortion, however, can be either static (permanent) or 
dynamic (fluctuating) 26 . Consider the case of a tetragonal distortion 
for which the Ca atoms in the axial positions move further away from 
the center of the octahedron; in that case, the schematic energy level 
diagram (Fig. 3) shows that the energy degeneracy of the t 2g set of 
orbitals is lifted and the resulting degenerate 3d xz and 3d yz orbitals 
then have the lowest energy. Thus, as these orbitals were filled, the 
electron -electron interactions are dominated by the equatorial Ca 
atoms located in the plane. To test this hypothesis, the electron 
density of a two-dimensional (2D) square net of Ca was calculated. 
Indeed, it was found that the electrons do occupy the Ca 3d xz and 3d yz 
orbitals and the net charge density is localized above and below the 
plane of the atoms (Fig. 3). As an extension to the tri-dimensional 
case 11 , the electron density maximum in the center of the cube can be 
obtained from the stacking of these 2D networks in a proper way. It is 
worth noting that the structure 27 and electron topology of tetragonal 
Cs IV 28 can be rationalized as the result of an s — > d transition 
followed by a static JT distortion, 

The JT theorem 26 does not provide guidance per se as whether the 
lattice distortion is either static or dynamic. To this end, constant- 
pressure/constant temperature (NPT) ab initio MD simulations 
using the pseudopotential planewave method 29 were carried out at 
several temperatures for the pressure of interest, namely, 42 GPa. A 
common observation arose from the calculations: the time-averaged 
structure remains cubic however significantly large lattice distortions 
are observed in the time domain. A representative trajectory 
obtained at 140 K is shown in Fig. 4: in this case, the 2X2X2 simu- 
lation cell fluctuates around the SC structure with an average cell 
length of 5.2 A. Several lattice distortion events indicated by arrows 
in Fig. 4 were captured in during the course of the simulation. Note 
that for these events, the cell axes and angles did not distorted ran- 
domly. In all cases, when one of the cell angles deviated from 90° to 
1 10°, a concomitant elongation of two cell axes to 5.4 A was recorded. 
It is conceivable that on time average, all cell axes and angles are 
sampled equally and the resulting cell maintains indeed the cubic 
symmetry. We note that, as a significant result, the lattice distorts 
with the approximate cell parameters oia = b = 5. 4 A, c = 5.1 A and a 
= P = 90°, y = 110° which, incidentally, are close to the cell para- 
meters of the predicted I^/amd structure 8 . The potential energy 
surface of the JT distortion is known to resemble a Mexican hat 25 , 
thus, inter- conversion dynamics as observed in the MD trajectories, 
sampled the potential energy surface without accessing the higher 
energy SC structure. The consequence of the dynamical fluctuations 
in dense Ca is that the cubic unit cell will present significant disorder 
in the [110] plane. As the atomic configuration is constantly chan- 
ging, it is considered as highly fluxional. In agreement with previous 
ab initio studies 9 " 12 the simple cubic structure is unstable at 0 K. 
Finite temperature lattice dynamics calculations, however, show 
the anharmonicity is reduced with increasing temperature and the 
simple cubic structure becomes dynamically stable at temperature 
higher than 200 K (Fig. 5). The observed small narrowing on the 
width of the (110) Bragg peak is consistent with the theoretical 
results. 

In summary, experimental and theoretical results show unam- 
biguously that the SC structure of Ca above 39 GPa is an "average" 
structure, presenting dynamic lattice distortions. There is no need to 
hypothesize that the SC structure should be energetically competitive 
to the predicted I^lamd structure at high temperatures 9 " 13 . In fact, 
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Figure 3 | A schematic representation of possible Jahn- Teller distortion of Ca in the SC structure under the effect of the crystal field (for explanation, 
see text). 



the energetically less favourable SC structure should not be sampled 
in the lattice fluctuations. At very low temperature from 39-42 GPa, 
due to structural fluctuations and local stress on the crystalline envir- 
onment, it is conceivable that a plethora of crystal structures of lower 
symmetry than the SC structure can be stabilized and observed 215 . 
The results presented although confirm the theoretical prediction 16 ' 30 
that the motions of Ca in phase III' are rather facile, however, the 
nature of structural distortion is different. In addition, an intuitive 



chemical principle was provided to explain on the origin of the large 
anharmonicity 30 . This insight should be applicable to the interpreta- 
tion of the structure and dynamics of elemental solids under high 
pressure. 

Methods 

Electronic calculations of the charge densities were performed with the VASP 
code 31 ' 32 using the projected augmented potential 33 for Ca with semi-core 3s and 3p 
orbitals treated as valence states. Constant pressure- constant volume molecular 
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Figure 4 | Temporal variation of the simulated cell angle (a) and lattice parameter (b). The changes in the cell dimensions can be interpreted as dynamic 
fluctuations of a cubic unit cell in the [110] plane, (c), (d) and (e) illustrates the "fluxional" distortions in the ob plane of the simple cubic cell. In a simple 
cubic cell, all the three crystal axes are formally equivalent. In the molecular dynamics simulation, one can labelled the axes as x, y, and z (4c). The figures 
show, schematically, that at a given instance, pair of the crystal axes, i.e. either the (x, y) or (x,z) axes are distorted simultaneously as expected for 
dynamical Jahn-Teller distortion of a cubic system. 



SCIENTIFIC REPORTS | 2 : 372 | DOI: 1 0.1 038/srep00372 



4 




dynamics calculations were performed with the Quantum Espresso code 
(www.pwscf.org) 34 using the pseudopotential reported in ref. 24 in which the semi- 
core electrons 3s and 3p together with those in the 4s orbital were treated as valence in 
the construction of the pseudopotential. A 2 X 2 X 2 supercell model was employed for 
the calculation. The Brillouin zone was sampled with a 4X4X4 /c-point mesh. 

Phonon dispersion curves were calculated at 0, 100, 150, and 200 K using the self- 
consistent ab initio lattice dynamics method 35 . Anharmonicity in the force constants 
are considered explicitly in the finite temperature simulation. The phonon calcula- 
tions employed a 4X4X4 supercell with a 4X4X4 k-point mesh. 
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